ON THE ENERGY OF INVISCID SINGULAR FLOWS 
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Abstract. It is known that the energy of a weak solution to the Eu- 
ler equation is conserved if it is slightly more regular than the Besov 

1/3 

space B 3 oa . When the singular set of the solution is (or belongs to) a 
smooth manifold, we derive various L p -space regularity criteria dimen- 
sionally equivalent to the critical one. In particular, if the singular set 
is a hypersurface the energy of u is conserved provided the one sided 
non-tangential limits to the surface exist and the non-tangential maximal 
function is L 3 integrable, while the maximal function of the pressure is 
L 3 / 2 integrable. The results directly apply to prove energy conservation 
of the classical vortex sheets in both 2D and 3D at least in those cases 
where the energy is finite. 



1. Introduction 

In this paper we study weak solutions to the Euler equations modeling 
evolution of inviscid fluid flows 

du 

(1) — + (u ■ V)u = -Vp, 

(2) V • u = 0. 

Here u is a divergence-free velocity field, and p is the internal pressure. The 
classical law of energy conservation 

\u(t)\ 2 dx = J \u \ 2 dx 

for smooth rapidly decaying solutions of (OQ) and © is an easy consequence 
of the antisymmetry of the nonlinear term. Weak solutions to CQ) are be- 
lieved to describe turbulent phenomena at large Renolds number in the in- 
ertial range of frequencies. The Kolmogorov-Obukhov power laws predict 
solutions to be |-Holder continuous in a statistically averaged sense. More- 
over, since the energy is not lost within the inertial range the energy flux 
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through inertial scales is to be proportional to the mean energy dissipation 
rate e ( IfTTlO . Experiments show that e is essentially independent of the vis- 
cosity coefficient. So, if in the limit of infinite Reynolds number turbulent 
solutions converge in some sense to weak solutions of the Euler solutions, 
then such solutions are expected to be on average energy dissipative. 

Onsager |[l4l stated that all (| + 5) -regular solutions conserve energy, 
and there may exist solutions exactly ^-regular that do not. The results of 
Eyink [9| followed by the work of Constantin, E and Titi flU give Onsager's 
hypothesis rigorous proof in the spaces Bl^ 5 , which measure Holder con- 
tinuity in the L 3 -space. An example of a vector field exhibited in (Hill sug- 
gests that the exponent | may indeed be critical, however no rigorous proof 
of this fact exists at the moment. An improvement upon by Duchon 
and Robert Q showed that some solutions conserve energy even in the 
Onsager-critical case. In recent paper [0 the criterion was established in 

1 /3 

the dimensionally optimal regularity class LfB 3 ' where cq signifies the 
decay 2 q \\ A ? w||| — > of the ^-derivatives of the dyadic parts. 

This present paper is motivated by the work of Caflisch, Klapper and 
Steele Q, where the authors obtain bi-Holder sufficient conditions for solu- 
tions with singularity set located on a smooth submanifold of IR n . Although 
these conditions are subcritical, they are more practical in applications, for 
example, to multifractal models of turbulence (see B21QI1)- However, other 
important classes of singular weak solutions such as vortex sheets remain 
unattainable by the results of [|2l [31 H). Indeed, classical analytic vortex 
sheets in 2D or in 3D fall exactly into the critical class Bl'^\Bl^. 

In this paper we study the energy law for solutions which exhibit or- 
ganized singular sets. Examples of singular set organization include time 
dependent families of submanifolds of W 1 and their locally finite unions. 
We obtain Onsager-critical criteria near such sets in terms of L p -spaces, 
which do not involve calculation of spacial Holder exponents. For instance, 
in the case of a 3D solution with point singularity s(t) at time t G [0, T] 
and s G C 3 / 5 ([0,T]) the energy of u is conserved provided u G L\Js!j 2 
near the curve s (see also application to viscous flows in |fT8l ). In higher 
dimension we use mixed L p -spaces relative to the singular manifold (see 
Theorem 13.21 and Section 13.11) . The case of hypersurface S(t) is treated 
separately in Section HI We will introduce the notion of a slit suitable for 
subsequent analysis. We assume that the velocity and pressure fields have 
non-tangential or normal limits and that the non-tangential maximal func- 
tions are integrable on the surface. As a consequence of weak formulation 
of the Euler equations, we show that all slits necessarily satisfy the kine- 
matic condition similar to that of a free surface, so that particles that are 
initially on the surface stay on the surface at all time (see Lemma |43T ). This 
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case is radically different from the lower dimensional case where no par- 
ticular evolution law is imposed by the equation. Our analysis shows that 
the energy of a solution u with a slit type of singularity is conserved pro- 
vided the non-tangential maximal functions of u and the pressure p belong 
to L 3 (S) and L 3 ^ 2 (S), respectively (see Theorem 14.41) . These conditions 
are verified for the classical 2D and 3D vortex sheets in Section [5] implying 
their energy conservation (under zero total circulation in 2D). 

Energy non-conservative weak solutions without any apparently orga- 
nized space singularities have long been constructed by Scheffer [fT6ll and 
Shnirelman IfTTl , and more recently by De Lellis and Szekelyhidi in ||5). 
Those belong to L 2 L 2 X and L^L 2 X , respectively, and therefore are consider- 
ably Onsager-supercritical. As we mentioned earlier the vector field con- 
sidered by Eyink [0 with non-vanishing energy flux belongs exactly to 
Bl^\Bl^ o . However, no weak solution with this initial condition is known 
to exist. The example serves to show that the traditional mollification ar- 
gument used to prove energy conservation is sharp. Again, one can show 
that it has no organized singularities. It is in fact locally nowhere in the 
energy-regular class B l J^ Q . 

Although we chose to use W 1 as a model case, the local nature of the 
arguments presented below allows us to apply the results to other boundary 
problems, such as periodic in all or some spacial directions. This will be 
especially useful in application to vortex sheets. 



2. Weak solutions and Regular sets 

Definition 2.1. A vector field u e C w ([0, T]; L 2 (R n )), (the space of weakly 
continuous functions), is a weak solution of the Euler equations with initial 
data u e L 2 (W n ) if for every ip G Cg°([0, T] x R n ) with V x • tp = and 
< t < T, we have 

(3) / u-ip— u -ip— / u-d s ijj= / / (u®u):Vip, 

JR n x{t} JR n x{0} JO JR n JO Jm. n 

and Vz • u(t) = in the sense of distributions. We define the operation : by 

A: B = Tr[AB]. 
It will be convenient to work with the associated pressure defined by 

n 

(4) p = - ^2 RiRk(uiu k ), 

l,k=l 

where Ri are the classical Riesz projections. With the use of p we can 
alternatively restate the definition of a weak solution without requiring V x • 
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ip = 0. Namely, 

(5) / u ■ ijj — u • ip — I I u ■ d s ijj 

JR n x{t} Jr«x{0} Jo Jm n 

/ (u (g) u) : Vijj + / / pdivip, 
'o Jn n Jo Jm™ 

holds for all ip E Cq°([0, T] x IR n ). Since the pressure is only a distribution, 
the pairing between p and div ip is to be understood accordingly. 

Based on the results of we introduce the global regularity class 1Z(M. n x 
/) consisting of vector fields u E L 3 (IR n x /) on a time interval / C [0, T] 
such that 

(6) lim -. — r / \u(x — y, t) — u(x, t)\ 3 dxdt = 0. 

y-+° \y\ Jr^xi 

For an open set U C IR n is an open set, we define 1Z(U x /) as the class of 
fields u such that u<p E 1Z{R n x J) for all <p E C%°(U), 

Alternatively, we could define 7KW 1 x /) using Littlewood-Paley decom- 
position over dyadic shells in the frequency space (see [fT9ll ) 



U = > L\ q U. 



J2 A 

q=0 

Thus, condition © is equivalent to 

(7) lim / 2 q \\A q u{t)\\ldt = 0. 

In this form the regularity class was introduced in and the energy con- 
servation was established. A similar but less time-optimal class was con- 
sidered in [7] as a direct improvement upon [4J. We remark that condition 
u E £ 3 ([0, T); Bl^) implies ©, where c stands to indicate 

lim 2 9 ||A (? w(t)|| 3 ; = 0. 

q— >oo 

Definition 2.2. Let u be a weak solution to the Euler equations. A point 
(xo, to) el"x [0, T] is called regular if there exists an open neighborhood 
U C K n of Xq and a relatively open interval I C [0, T] containing t such 
that m G 7£(?7 x /). An open set D C IR n x [0, T] is regular if every point 
in it is regular. The set S of all irregular points is called the singular set of 
u. 

The main purpose of this section is to prove the following local energy 
balance relation inside every regular set. For a set A C M" x [0, T] we 
denote by A(t) the slice AnPx {£}. 
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Lemma 2.3. Let D be a regular set of a weak solution u. Then for every 
4> G Cq°(D) one has 

(8) / \u\ 2 <p- i |w| 2 0- / \u\ 2 d t <p= \ (\u\ 2 + 2p)u- V0, 

JD(t") Jd(V) JD JD 

for alii!, t" G [0,T]. 

Before we prove this lemma, we need to take another seemingly obvious 
but not entirely straightforward step by showing that one can substitute a 
mollified in space solution u into © as a test function. This fact is not 
so straightforward since a priori u may not have sufficient time regularity. 
The difficulty has been removed in a similar situation in [fi3l by considering 
mollification both in space and time, however in our case such mollification 
would introduce unnecessary technical obscurity. So, let us fix a mollifier 
h G C™(M. n ) with / h = 1 and h = outside the unit ball. Denote 

us(x,t) = / hs(y)u(x-y,t)dy, 
h s (y)=5- n h(y5- 1 ). 
We prove the following lemma. 

Lemma 2.4. Let u be a weak solution. Then for each fixed 5 > 0, u§ '■ 
[0, T] — > W s,q is absolutely continuous for all s > and q > 2, and 
moreover 



(9) d t u s = -V • (u ® u) s - V Ps , 

fora.e. t 6 [0,T]. 

Proof. Substituting test-functions of the form 

where f3 G C °°(0, T) and ip G L7 °°(M n ) into © we obtain 

dt / u(t) ■ i) = / u(t) <g> u(t) :Vip+ p(t) divip 



in the distributional sense. Hence, since u is weakly continuous 

j u(t)-ip = J u(0)-ip + J J u(s)®u(s) :Vtl)ds + J J p(s)divipds } 
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for all < t < T. Let if>i denote the coordinate components of if). Taking 
the Fourier transform we obtain 



u(t) I u(0) -ijj-i I I (u lUk ) A (Z, s^MOdCds 

-iff p(Z,s)ZMOd£ds, 



assuming the usual summation convention. Let us notice that (uiUk) A and p 
are continuous and bounded functions of £ 7^ for every s. Let E t denote 
the common Lebesgue set of u(t) and u(0) not containing the origin, so that 
|IR n \S t I = 0. Denote by &,-(£), j = 1, n the vectors of the standard unit 
basis. For every j and £ G E t we apply the previous identity to a sequence 

of functions if) such that if) n (£) — > e 3 -(£)5o(- — £), where <5 is the Dirac mass. 
We obtain 

Uj(C,t)-e a (0 =ttj(£,0)-e a (f)-i / {ujU h ) A (s, / p(£, s)&cis, 

Jo Jo 

for all i G [0, T] and £ G E t . Thus, the identity 

(10) u(t) = n - / [V • (« ® u) + Vp]ds 

holds in the sense of distributions for all t G [0, T]. Mollifying (PTOl) with foj 
we obtain 

(11) u s {t)=u s {0)- f [V • {u®u) s + Vps}ds 



Jo 

for all t G [0, T]. Since u(0) G L 2 we have u 5 (0) G l^ 5 ' 9 for all s > and 
q > 2, and since w ® u G LfL l x , we have V ■ (it <g> u)^ Vps G LfW*> q , 
for all s > and q > 1. This proves the lemma. □ 

Let us denote 

£ c =(-c,c) n 

Proof of Lemma \2J\ First let us observe that p G Lf^(D). Indeed, for a 
compact subset if C D let e > be such that K + B e C D. Let a G 
Cg°(D) be such that a = 1 on F + £ e/2 and a = on D\K + B e . Then 

(12) p = RiRj(uiUjip) + RiRj(uiUj(l — a)). 

Since UiUja G L 3 / 2 , so is the first term in (fT2l) . The second term belongs to 
L°°(K) since 1 — « = on K + -E> e /2> UiUj G L 1 and the kernel of RiRj is 
bounded away from the e/2-neighborhood of the origin. This observation 
justifies the pressure integral in ®. 
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Using partition of unity over the support of <p we reduce the lemma to the 
case D = U x I, where U is an open ball. So, suppose 4> £ C^(U x /). 
Choose 5 > so small that 



supp(0(-,t) + B So ) C K C U, 
for all t E I. Let us now use © with 5 < 5 . We obtain 

d t u s ■ u S (j> = ^d t (\u 5 \ 2 (f)) - -\u s \ 2 d t (f). 
Integrating in time on [t', t"} C [0, T] we obtain 

(13) / M 2 0- / \u s \ 2 (f)- / \u s \ 2 d t t 

JUx{t"} JUx{t'} JUx[t',t"] 



= 2 (u®u) 5 : V(u S (f)) + 2 p 5 u s ■ V0. 

Jux[t',t"] Ju-x.lt', t"\ 

Notice that the time integration is in fact happening on the interval [t', t"] Hi. 
So, we can pass to the limit as 5 — > on the left hand side and in the pressure 
term. The nonlinear term will be treated similar to [4J. First, consider a 
scalar (5 E C^°(U) with /3 = 1 on K. We can then replace u by u(3 under 
the integrals of (fl"3l) . Without further change of notation we simply assume 
u E TZ(R n x I). We have 

(u®u) s : V(m 5 0) = 2 / (u <g> u) s : (u s <g> V0) 

!7x[t',t"] Jux[t',t"} 



Clearly, we can pass to the limit 



+ 2 (u®u) s : V{u s )(j). 

>Ux[t',t"] 



(14) 2/ (m®m) 5 : (u s ® V0) ->• 2 / • V0. 

JUx[t',t"] JUx[t',t"] 

Let us observe the following identity 

(15) (u ® u)s = rs(u, u) — (u — us) ®(u — us) + us <g> -u<5, 
where 

r«5(w,«)(a;,t) = / h s (y)(u(x-y,t)-u(x,t))®(u(x-y,t)-u(x,t))dy. 
Notice 

u(:r) - u s (x) = / h s (y)(u(x) - u(x - y))dy, 



and 

V«j(i) = - / (Vh) s (y) <g> - m(x - 
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So, we can estimate using Holder and Minkowski inequalities 



(u - us) ® {u — u 5 ) : V(u s )(j) 

< I / hs(y) / \u(x, t) — u(x — y, t)\ 3 dxdtdy J 

x (f l\(^ h )s(y)\ [ \u(x,t) - u(x - y,t)\ 3 dxdtdy] 

~w(J ^ hs ^ dy ) m (/ I^IKV^^I^V 73 < o(5) -> 0. 
Similarly, the term with vanishes as well. Finally, 

2 J us ® us : V(us)(f) = — J \us\ 2 us ■ V(f> ^ ~ J \ u \ 2u ' ^<P- 
This adds up with (fl4l) to produce the corresponding term in ([8]). □ 

3. Low-dimensional singular sets 

Definition 3.1. We say that a set S C W 1 x [0, T] admits a k-dimensional 
C 1,l -cover if for every point (xo,t ) in the space-time there is an open 
neighborhood U of x in M n and a relatively open subinterval J C [0,T] 
containing t for which there exists a family of C 1 -diffeomorphisms 

(16) <f t :U^B h tel, 

satisfying the following conditions 

(a) S(t) n U c ¥>r x (R* x {0} n - fe n Bi), for all t G I; 

(b) There is C > such that 

sup \tp t >(x) - ft"(x)\ < C\t' - t"| 7 , 
xeu 

for all f , t" G /; 

(c) sup^g^gj |V x p t (a;)| < C. 

Theorem 3.2. Le? u G L 3 (IR n x [0,T]) &e a weafc solution to the Euler 
equation on the time interval [0, T]. Then u conserves energy provided the 
singular set S of u admits a k-dimensional C^^-cover and u G L?Ll oc , 
where the values ofj,n,k,q > satisfy 

(IV) 7 > t £ tt, n > k + 2, q > 3- ?/ ~ " 



(q — 2) (n — k) ' ' n — — 1 
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Proof. We claim that in order to prove Theorem 13 .21 it suffices to show that 
for every coordinate chart U x I and scalar test-function <\> E C^°(U) inde- 
pendent of time one has the following identity 

(18) / \u\ 2 (p- I \u\ 2 (p= I (\u\ 2 + 2p)u- V0, 

JUx{t"} JUx{t'} JUx[t',t"] 

for all t" E I. Indeed, if this is the case, we fix an arbitrary smooth <p 
with supp(0) C B R , and t E [0, T]. By compactness we can find a finite 
collection of charts UiXli,i = 1, M so that all Ij's contain t and Ui's cover 
B R . Put I = nf^/j. Consider a partition of unity {gt}^ subordinate to 
the cover, so that supp gi C Ui and J2i=i 9i = 1 on Br. Since we have (fl"8l 
for any and t', t" E Iq summing up over i we obtain (fT8l) for the given 
itself. The above construction is carried out for every t E [0, T\. Thus, we 
can find a finite cover of [0, T] by intervals such as I , and as a consequence 
obtain C[8]) for all t', t" E [0, T}. Letting = 4> (x/R), where O = 1 on 
£?i and O — on B 2 , and letting R — > oo we see that the right hand side of 
(fl"8~l) vanishes and we arrive at the desired energy equality. 

We will prove (fl"8l with the use of Lemma 12.3 L but first we need to in- 
troduce a cut-off of the singular sets S(t) n U. Let (ft : U — > i?i be the 
coordinate map, for t E I. Denote / = [a, 6]. If <^ t is not defined at a or 
b, then £ is not that point. In this case we can consider a slightly shorter 
interval / still containing t and so that ip t is defined at both ends. Let us 
define an extension of <p t as follows 

{ip a , t <a 
<pt, a<t <b 
(p b , t>b 

Notice that (p t still satisfies condition (b) of Definition 13. H on the entire real 
line. Let /3(r) be a mollifier. Define 

Jr 

Let us notice the following approximation inequalities: 

(20) sup \(ft, e ( x ) ~ V 2 *^)! < Ce 1 ] 
xeu,tei 

(21) sup \d m>e (x)\ <Ge<- x ; 

x£U,t£l 

(22) sup \V x <pt, e (x)\<C. 

xGU,tei 
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Let us fix a non-negative function 77 E C^°(W l ~ k ) with 77 = 1 on B 2 c and 
77 = on M n_A: \i?3 C '. We consider the following cut-off function 

( 1 1 

Xe(t,X) = 1-77 l—(ip t)€ (x)) k +l, —((pt,e(x))n 

for x E U and t E I. Notice that as e — > x e — > 1 for all t and a.e. x. 
Furthermore, due to (l20l) . supp(x £ 0) does not intersect the set S on the time 
interval /. Finally, put 

<Pe = Xe<P- 

Due to regularity of u away from S, Lemma [2731 applied to produce 
(23) / |«| 2 e - / M 2 e - / \u\ 2 d t <p e 

JUx{t"} JUx{t'} JUx[t',t"] 



[ (\u\ 2 + 2p)u- V0 £ . 

JuxU't"} 



'Ux[t',t"] 

Let us examine the terms in the limit as e — > 0. Clearly, the first two terms 
on the right hand side will converge to their natural limits. As to the third 
term, we have d t (f) e = (f)d t Xe, and 

dtXe = -— dtiwA^Av (—(¥tA x ))k+i,---,—(vtA x ))n) ■ 

j=k+l ^ ' 

Notice that d t Xe is supported on the set 

U n ((p tte )- 1 (R k x [-3Ce 7 , 3Ce 7 ]^ fe n B x ), 
which is a subset of 

A e = U n (y^rW x [-4Ce 7 , 4Ce 7 ] n - fc n Bi). 
We have |A e | ~ e (n " fc)7 . In view of <(T3b and (ED we obtain 



(24) 



\u\ 2 4>d t x ( 



Ux[t',t") 



\ A \{q-2)/q ft" / f \ 2/9 

< i-^ / / dt 



t" / p \ 2/9 



(25) '^0. 
Let us now examine the right hand side of (|23l) . We have 

Clearly we can pass to the limit in the integral containing the second term. 
As to the first term we have 

1 n ( 1 1 \ 

VXe = Yl V &tA X ))AV {—( l fit,e( X ))k+l,---,—( l ft,e( X ))nj , 

j=k+l ^ ' 
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which is supported on the set A e . Thus, 



(26) 



/ (\u\ 2 + 2p)u- V0 e 

Jux[t',t"] 



\A |(<?- 3 )/9 /■*" ft \ 3/q 

< / ( / \u\ q dx \ dt 



A, 



rt" f r \3/g 

(27) ~ / / \u\ q dx dt — > 0. 



This finishes the proof of Theorem 13 .21 □ 

3.1. An Onsager-critical improvement. Let us consider physical units of 
velocity - U, length - L and time - T. Then the dimension of the regu- 
larity space 1Z is T^f/L 2 ^ - . We call functional spaces of this dimension 
Onsager-critical. In the case of point singularities, i.e. k = 0, Theorem 13 .21 

q 3n 

yields the Onsager-critical condition u G L^L™- 1 with 7 being at least 
^pg. Under these circumstances we expect our result to be optimal. How- 

3(n-fc) 

ever, this is not the case if k > 0, since the dimension of L^L™-*- 1 is 

n(ri — fc — 1) 

T l ' 3 UL 3(n-k) _ Onsager-critical spaces for > can be defined using 
mixed L p spaces relative to the slices S(t). Assuming that each S(t) is a k- 
dimensional smooth submanifold of R™ we consider local normal fiber bun- 
dle 5 ,± (t). Thus, each fiber S- L (x, t) is a 7-smooth in time local tile orthog- 
onal to the surface S(t). We can now define the local space u E L^L P S L 9 S± 
by requiring over coordinate neighborhood U x I the condition 

[([ ([ \u(x,y,t)\«dar k (y)Y 9 da*(x)\ dt < 00, 

Ji \Js(t)nU KJs^ix^nu J ) 

where da t indicates the surface measure of the corresponding dimension. 

3(n-fc) 

Notice that the space L 3 L|L^ fc_1 is in fact Onsager-critical. In general, 
Theorem 13 .21 can be restated by requiring 

(28) u G L 3 t L 3 s L g s± 

under the same assumptions on n, k, 7, q. In particular, we obtain energy 
conservation if 

3(n-fc) 3 

(29) u G L 3 L|L^ fe_1 and 7 > 



n - k + 2 

In order to reprove Theorem 13.21 under new condition (|28l) one has to 
simply apply the Holder inequality in (|24|) and in (|26l) only to the integrals 
over S ± (x 1 1), the rest of the argument being the same. We leave details for 
the reader. 



12 



ROMAN SHVYDKOY 



3.2. Other extensions. Since our argument is local, it is readily extendable 
to the case of locally finite union of singular sets. Specifically, suppose that 
in every coordinate neighborhood V = U x I 

N v 

(30) S = (J Sj, 

i=i 

where S/s are /cj-dimensionally C 7j ^-covered in V. We can use the prod- 
uct of cut-offs 

N v 

Xe = Xl 
3=1 

to run the argument. The conclusions of Theorem 13.21 remains true under 
the corresponding assumptions on u locally near each Sj. The result of 
Section [3TT1 can be modified similarly. 

4. The case of hypersurface: slits 

In this section we will study the case k = n — 1. We will assume special 
geometric properties of the singular set S. Namely, let S be a C l -family 
of closed orientable C 2 -submanifolds of R™. For every (x , t ) E S there 
exist U, I and a local parametrization r = r(y, t) of S(t) D U for all t E I, 
where r E C|' t , and y = (y a , y n _i) G Let z7(x, t) be the positively 

oriented unit normal to S(t). We consider a coordinate system on a smaller 
neighborhood that is most suitable for dealing with normal limits. For e > 
small we define 

A(v, Vn) = r(y, t) + e y n v(r(y, t),t), 

for \y n \ < 1. Since S is sufficiently smooth, this defines a diffeomorphism 
of B™ onto an open neighborhood U(t) with S(t) n C/(t) = S(t) n C/ for all 
t E I. It will be convenient in the future to deal with [/ independent of t. 
So, reducing the time interval if necessary we can find a new neighborhood 
U CU(t) for all t E I, such that 

I)"" 1 x cO) C [/ C ^((-1, l)"" 1 x (-e 2 , e 2 )) 

for alH G / and some c 2 > c\ > 0. The direct product V — U x I along 
with the map ip t = ipt define a new coordinate chart containing (xo, t ). 
Let us also define the normal segments for every (x, t) E V: 

V+(x,t) = (x,t) + z7(x,t)[0,ei] 
r_(x,t) = (x,t)+u(x,t)[-e 1 ,0}. 
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We may further truncate the segments to ensure that for some open neigh- 
borhood W of S we have \J S T± C W . For a function or field / on W we 
denote by f± : S — > K the normal maximal function defined by 



f±(x,t)= sup |/(a/,t)|, 

x'6T±(a!,t) 



and by / ± the limits 



/±(x,t) = lim J(x',t), 

x'— »x,x'£r± (x,t) 

if the latter exist. 

We now introduce a measure on each S(t) whose role will be clear in a 
moment. We start by defining it locally on every chart U H S(t). For this 
purpose let us fix a scalar-valued function H(x, t) G C 1 with level surface 
{H(x,t) = 0} = S(t) n E7 for all £ G /, and such that V X H ^ agrees 
with v. For instance, H = ((ft(x)) n . Let us consider the measure 

(31) dtf(x) = -^Ld* t (x) 

| v x H I 

where da t {x) is the surface measure of S(t). Notice the following identities 

H(r(y,t),t) = Q 

d t H{x,t) = -d t r(y,t)-V x H(x,t), 
where x = r(y,t) G S(t) fl [/. Thus, in local coordinates, 

dnY(x) = -d t r(y,t) ■ u(r(y,t),t)J t (y)dy, 

where Jt(y) is the volume element. We see that the definition of dfj^ is 
independent of H . Yet (|3TT) shows that it is also independent of particular 
parametrization of S(t). Now, let / G C (S(t)) be a continuous function 
with compact support on S(t). Arguing as in the proof of Theorem 13 .21 we 
find a finite cover of supp(/) by {f/, }*^ with the corresponding partition 
of unity {gi}^ x over supp(/). Define 

(32) / fdjjL t {x) = fgid/if'ix). 
Js(t) i=1 Js(t) 

This is a well-defined measure over S(t). For instance, if S(t) is given by 
the graph of a periodic in spacial variables function x n = z(x X: ...,x n ^x,t), 
then 

dfi t = -d t z(xx, ...,x n -i,t)dx 1 ...dx n - 1 . 
The measure dfj, t arises naturally in the following calculation. Let us fix 
a coordinate chart (V, <p t ) a s above, define rj as in the previous section with 
k = n — 1, and denote 

(33) X e(x,t) = l-r ] (e-\ Vt (x)) n ), 



14 ROMAN SHVYDKOY 

for e < e\. 

Lemma 4.1. Let f : V — > R and u : V ^ W 1 be such that the limits 
f±(x, t) and u±(x, t) exist for a.e. t E I and a.e. x G S(t) with respect to 
da t , and f±, u* ± G L 1 (da t dt). Then 



lim / fdtXedxdt 



US(t) 



(34) 
and 

(35) lim / u ■ V ' x Xedxdt = (w+ — W-) • Vdo T dr. 

Jvo ' .// ./s(t) 

iVoo/ Let us denote H(x, t) = ((p t (x)) n . To prove (l34l) let us observe 

/ fd t Xedxdt = - [ e- 1 r ] '(e- 1 H(x,t))d t H(x,t)f(x,t)dxdt. 
Jv Jv 

As a guiding point we recall the classical microlocal limit 
- / gdx^ -TWT7] da - 

e J0<H<e JH=0 I V-H I 

By changing the variables we obtain the integral 



fdtXedxdt 



F e (y,t)dydt, 



where 

Fe(y,t) 

and 



fMy,yn),t)d t H(ip t (y,y n ),t)e 77 (e y n ) 



\Vn\<fil 



fLt{y,y n )dy n dt, 



det 



DA 



Dy 



My) 

Given our choice of H we have 

(36) H(ijj t (y,y n ),t) = y n , 

(37) d t H + V x H-dMy,y n ) = 0. 

So, as y n — ► we obtain 

d t H^-d t r(y,t)-V x H(r(y,t)) 
uniformly in y G 5™ . Moreover, 

Mv) e o 
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Using that V x Ii(x ) t) = e 1 z7(x, t) we obtain the uniform convergence 

d t HVt t -> -d t r ■ uJ t (y). 

Let us observe now that as e gets sufficiently small, we have 4>t(y, y n ) £ 
F sg n(y n ){ipt{y, 0),t) for all y E B^ 1 , and ip t (y, y n ) approaches the surface 
orthogonally. The condition f± E L l (da t dt) implies that all F € have a 
common integrable majorant. This enables us to pass to the limit and arrive 
at (|34|> . The proof of (03]> is similar. □ 

Definition 4.2. Let u be a weak solution to the Euler equations. The surface 
S is called a slit of k if 

1) The limits u±,p± exist for a.e. t E [0, T] and a.e. x E S(t), 

2) u*± E L 2 (da t dt)\ oc andp± E L 1 (dcr t dt)i oc 

Lemma 4.3. Le? u be a weak solution to the Euler equations, and S be a 
slit. Then the following is true: 

1) u + • V = u_ ■ V := u v and p + = p_ for a.e. t E [0, T] and a.e. 
x E S(t); 

2) dfi t + u v da t = Ofor a.e. t on the set u + ^ m_. 

Proof. As before we reduce the statements of the lemma to the local coordi- 
nate neighborhood V = U x I defined earlier. Let us consider an arbitrary 
scalar function g E Cq°(V). From the divergence-free condition on u we 
obtain 



u ■ V(gxe) = 0. 

v 

Letting e->0we obtain from Lemma |4~T1 

/ u ■ Vg +jj g( u + ■ v — ■ u)da T dr = 0. 
Jv J I Js{t) 

Using the divergence-free condition again and the free choice of g we obtain 

(38) u + ■ V = u- ■ V. 

Consider an arbitrary vector-valued function a E Cq(V), and ip = ax € - 
By continuity, the regularity of ip is sufficient to substitute ip into ©. We 
obtain the following identity: 

— / u-d T aXe— / u-ad T Xe= \ (u <8> u) : Va Xe + / {u ■ a){u ■ Vxe) 
Jv Jv Jv Jv 

+ \ p{a- Vxe + Xediva). 
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Using (1341) and (1351) we obtain in the limit as e — > 



— u ■ d T a — I I {u+ — u~) ■ adp T dr 

>V JlJs(r) 



u®u): Va + / / ( M + — w_) • au v da T dr 

v J I J S(r) 



+ 11 (p+ — P-)a u da T dT + / pdiva. 

J I JSM JV 



'S(.r) 

Using the identity for the weak solutions © with ip = a we see that only 
the boundary terms remain: 



(u + — uJ) ■ adij, T dT =11 (w+ — u -) ' du u da T dr 
lJS(r) JiJs(t) 



(p+ — P-)a u dcr T dT. 

US(r) 

Let us choose a of the form a = ug, where g G Cq(V). Using (|38|) we have 

(p+ — P-)gda T dr = 0. 



7 JS(r) 

This readily implies p + = p_ a.e. Going back to the previous identity we 
notice that 2) holds as well due to arbitrariness of g. □ 

Theorem 4.4. Suppose that u G L 3 (M n x [0, T]) is a weak solution to the 
Euler equations and the singular set Sofu is a slit. Suppose further that 
u± G L 3 (da t dt)i oc , p* ± G L 3 ^ 2 (da t dt)i oc . Then u conserves energy. 

In view of our discussion in Section [3T| we notice that the conditions of 
Theorem 14.41 are Onsager-critical. We therefore expect these conditions to 
be optimal as far as our argument in concerned. 

Proof. As in the proof of Theorem l3.2l we reduce the problem to proving the 
local energy equality (TT8T) . As before H (x, t) = (ip t (x)) n and <p e is defined 
by (1331) . The regularity of u away from the slit S enables us to use Lemma 
12.31 with <j> e . Using the results of Lemmas 14.11 and 14.31 we can pass to the 
limit as e — > and obtain 

u\ 2 d t <f) e -> / \u\ 2 <j) 

Ux[t',t"} JUx[t',t"] 
t" r 

2 |„, 1 2\ 



+ / / (\u+\ — \u~\ )4>dfi T dT, 

If JS(t) 
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and 



/ (\u\ 2 + 2p)u- V(f) e -> / (|w| 2 + 2p)n- V0 

JUx[t',t"] JUx[t',t"] 

+ / / (l M +| 2 ~~ \u-\ 2 )u v da T dT 
Jt' Js(t) 
,-t 



2 / / (p+ — P-)u u da T dr. 
Jt' JsM 



'S(r) 

According to Lemma l43l the surface integral terms sum up to zero, and (fl"8l) 
follows. □ 

Arguing as in Section l3?2~1 we can include the result of Theorem l4.4l in ob- 
taining more general singular set configurations. Thus, the union (|3Q|) may 
involve finitely many slits accompanied by the corresponding conditions on 
u and p. 

We remark that one can also state the conditions of Theorem 14 . 41 and Def- 
inition 14.21 in terms of more conventional non-tangential limits and max- 
imal functions. It would be interesting to know whether the condition 
u± e L 2q (S) automatically implies^ G L q (S). 

5. Energy of vortex sheets 

Naturally, the conditions of Theorem 14 .41 apply to vortex sheet solutions. 
Vortex sheets in the classical sense (as opposed to those defined by Delort 
O) are singular solutions to the Euler equations with vorticity concentrated 
on a hypersurface (see lfT5lO . For notational convenience we will consider 
the two dimensional case, although all what follows holds true in three di- 
mensions as well. In 2D a vortex sheet is described by the graph of a regu- 
lar function ((a, t) = (a, h(a, t)) and vorticity density 7 = 7(0, t) on the 
graph. Typically, one assumes 2n -periodicity on h and 7. Thus, in complex 
variable notation the velocity field off the sheet is given by the Biot-Savart 
law 

u ^ z - ,t > = J Cot I 2 J ' 1 ^ a ' t > da - 

Provided 7 has enough smoothness on a time interval [0, T], the standard 
potential theoretical considerations imply that u E L^L™, the non-tangential, 
and hence normal, limits exist and are given by 

u± (a, t) = — -pv r cot ( CIM-CK^A t s da , T , t , g 

4ttz J_ n V 2 / 
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where s is the unit tangent vector oriented in the positive direction of the x- 
axis. The pressure can be recovered from Bernoulli's function, and is given 
by the double-layer potential formula 

P=-\H 2 + \v(\u + \ 2 -\u_\ 2 ). 

From the classical jump relations for the double-layer potential V we con- 
clude that the limits p± exist, p + = p_ = |(|m + | 2 + |u_| 2 ) and p* ± 6 
L q (da t dt)\ oc for all 1 < q < oo. Thus, according to Definition 14.21 the clas- 
sical vortex sheet is a slit. The equation 2) in Lemma 1431 is nothing but the 
well-known evolution law of the sheet: 

dth = -U x d a h + U 2 , 

where U = + uJ). In order for the total kinetic energy of the vortex 
sheet to be finite we assume vanishing of the total circulation: 




7(a, t)da = 0. 



Under this condition, u £ L°°L 2 . By interpolation with u 6 L^L 00 we 
obtain u E L 3 L 3 . Therefore, the conditions of Theorem 14.41 are satisfied 
and we arrive at the following corollary. 

Corollary 5.1. Suppose that e C°°([0, T] x [— 7r,7r]), and the total 
circulation ofj is zero. Then the energy of the vortex sheet is conserved. 

Vortex sheets of this nature are known to exist in 2D and 3D locally in 
time in spaces of functions that admit analytic extension to a complex strip 
(see [HI |20]|). In general, the global existence is precluded by occurrence 
of the roll-up singularity (see lfT2lO . The conditions on Cauchy data stated 
in [|20l that guarantee local existence allow for sheets with zero circulation. 
Thus, Corollary 15.11 applies to a variety of existing vortex sheets. However, 
the proof of Theorem 14.41 applies to obtain local energy balance relation for 
sheets with infinite energy as well. 



References 

[1] Russel E. Caflisch. Long time existence and singularity formation for vortex sheets. 
In Vortex methods (Los Angeles, CA, 1987), volume 1360 of Lecture Notes in Math., 
pages 1-8. Springer, Berlin, 1988. 

[2] Russel E. Caflisch, Isaac Klapper, and Gregory Steele. Remarks on singularities, di- 
mension and energy dissipation for ideal hydrodynamics and MHD. Comm. Math. 
Phys., 184(2):443^155, 1997. 

[3] A. Cheskidov, P. Constantin, S. Friedlander, and R. Shvydkoy. Energy conservation 
and Onsager's conjecture for the Euler equations, (preprint). 



ON THE ENERGY OF INVISCID SINGULAR FLOWS 



19 



[4] Peter Constantin, Weinan E, and Edriss S. Titi. Onsager's conjecture on the energy 
conservation for solutions of Euler's equation. Comm. Math. Phys., 165(l):207-209, 
1994. 

[5] Camillo De Lellis and Szekelyhidi Laszlo. The Euler equations as a differential in- 
clusion, (preprint). 

[6] Jean-Marc Delort. Existence de nappes de tourbillon en dimension deux. J. Amer. 
Math. Soc, 4(3):553-586, 1991. 

[7] Jean Duchon and Raoul Robert. Inertial energy dissipation for weak solutions of in- 
compressible Euler and Navier-Stokes equations. Nonlinearity, 13(l):249-255, 2000. 

[8] Jens Eggers and Siegfried Grossmann. Origin of the Obukhov scaling relation in tur- 
bulence. Phys. Lett. A, 153(1):12-15, 1991. 

[9] Gregory L. Eyink. Energy dissipation without viscosity in ideal hydrodynamics. I. 
Fourier analysis and local energy transfer. Phys. D, 78(3-4):222-240, 1994. 
[10] U. Frisch and G. Parisi. On the singularity structure of fully developed turbulence. In 
M. Ghil, R. Benzi, and G. Parisi, editors, Turbulence and predictability in geophysical 
fluid dynamics and climate dynamics, Proc. International Summer School of Physics 
"Enrico Fermi", pages 84-87. Amsterdam: North-Holland, 1985. 
[11] Uriel Frisch. Turbulence. Cambridge University Press, Cambridge, 1995. The legacy 
of A. N. Kolmogorov. 

[12] Robert Krasny. A study of singularity formation in a vortex sheet by the point- vortex 

approximation. J. Fluid Mech., 167:65-93, 1986. 
[13] Igor Kukavica. Role of the pressure for validity of the energy equality for solutions of 

the Navier-Stokes equation. J. Dynam. Differential Equations, 18(2):46 1-482, 2006. 
[14] L. Onsager. Statistical hydrodynamics. Nuovo Cimento (9), 6(Supplemento, 2(Con- 

vegno Internazionale di Meccanica Statistica)):279-287, 1949. 
[15] P. G. Saffman. Vortex dynamics. Cambridge Monographs on Mechanics and Applied 

Mathematics. Cambridge University Press, New York, 1992. 
[16] Vladimir Scheffer. An inviscid flow with compact support in space-time. J. Geom. 

Anal., 3(4):343-401, 1993. 
[17] A. Shnirelman. On the nonuniqueness of weak solution of the Euler equation. Comm. 

PureAppl. Math., 50(12): 1261-1286, 1997. 
[18] R. Shvydkoy. A geometric condition implying energy equality for solutions of 3D 

Navier-Stokes equation, (preprint). 
[19] Elias M. Stein. Harmonic analysis: real-variable methods, orthogonality, and oscil- 
latory integrals, volume 43 of Princeton Mathematical Series. Princeton University 

Press, Princeton, NJ, 1993. With the assistance of Timothy S. Murphy, Monographs 

in Harmonic Analysis, III. 
[20] C. Sulem, P.-L. Sulem, C. Bardos, and U. Frisch. Finite time analyticity for the two- 

and three-dimensional Kelvin-Helmholtz instability. Comm. Math. Phys., 80(4):485- 

516, 1981. 



Department of Mathematics, Stat, and Comp. Sci., 851 S Morgan St., 
M/C 249, University of Illinois, Chicago, IL 60607 
E-mail address: shvydkoy@math . uic . edu 



